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ABSTRACT 

We compute masses and densities for ten periodic comets with known sizes: lP/Halley, 
2P/Encke, 6P/d'Arrest, 9P/Tempel 1, lOP/Tempel 2, 19P/Borrelly, 22P/Kopff, 
46P/Wirtanen, 67P/Churyumov-Gerasimcnko and 81P/Wild 2. The method follows 
the one developed by Packman and colleagues (Rickman 1986, 1989; Rickman et al. 
1987), which is based on the gas production curve and on the change in the orbital pe- 
riod due to the non-gravitational force. The gas production curve is inferred from the 
visual lightcurve. We found that the computed masses cover more than three orders 
of magnitude: ~ (0.3 - 400) x 10 12 kg. The computed densities are in all cases very 
low ( H 0.8 g cm -3 ), with an average value of 0.4 g cm -3 , in agreement with previous 
results and models of the cometary nucleus depicting it as a very porous object. The 
computed comet densities turn out to be the lowest among the different populations of 
solar system minor bodies, in particular as compared to those of near-Earth asteroids 
(NEAs). We conclude that the model applied in this work, in spite of its simplicity (as 
compared to more sophisticated thcrmophysical models applied to very few comets) , is 
useful for a statistical approach to the mean density of the cometary nuclei. However, 
we cannot assess from this simple model if there is a real dispersion among the bulk 
densities of comets that could tell us about differences in physical structure (porosity) 
and/or chemical composition. 
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1 INTRODUCTION 

Comets have a great cosmogonical interest since they are 
thought to be the relics of planet formation. They are ice- 
rich bodies formed beyond the snowline, i.e. the distance 
from the Sun at which water ice could condense. Comet ma- 
terial is presumably the most pristine so far known in the 
solar system that may keep the record of the grain accretion 
process in the protoplanetary disk. Micrometer-sized grains 
will tend to collect after gentle collisions into loose aggre- 
gates. Since comets are very small bodies, their self-gravity 
was never large enough to compact the material, so they 
are expected to conserve their primordial porous structure. 
Therefore, the comet nucleus is regarded as a very fragile 
and low-density structure. Donn (1963) estimated that the 
mixture of ices and meteoric matter would form an aggre- 
gate with a density of a few tenths g cm -3 . Later Weissman 
(1986) proposed a nuclear structure consisting of aggregates 
of various sizes that he called the 'primordial rubble-pile' 
model. 

The comet bulk density is a key parameter to learn 
about the geochemistry and the degree of compactness of 
the comet material. Masses and sizes of comets are nec- 
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essary to derive their bulk densities. The determination of 
these two physical parameters is not trivial. The main diffi- 
culty to determine the comet size is that comets are usually 
shrouded by a coma of dust and gas. Nevertheless, great ef- 
forts have been made to observe Jupiter family comets near 
aphelion where they are inactive or have only residual ac- 
tivity, so estimates of comet sizes have been possible for a 
large sample of Jupiter family comets (e.g. Fernandez et al. 
1999, Tancredi et al. 2000, 2006). It has also been possi- 
ble to obtain very accurate sizes of a few comets from close 
up images taken from spacecraft during flyby encounters. 
These are the cases of lP/Halley (e.g. Keller et al. 1987), 
19P/Borrelly (e.g. Soderblom et al. 2002), 81P/Wild 2 (e.g. 
Brownlee et al. 2004, Duxbury et al. 2004, Howington-Kraus 
et al. 2005), and 9P/Tempel 1 (e.g. A'Hearn et al. 2005). 

The determination of comet masses has been an even 
more difficult task. In 1805 Laplace could set an upper limit 
of 1/5000 of the Earth mass for the mass of comet Lexell by 
noting that the Earth's motion did not suffer any perceptible 
perturbation after a close encounter with the comet in 1770. 
From then new attempts to derive comet masses led to even 
smaller masses showing that comets were indeed very small 
bodies. Vorontsov-Velyaminov (1946) derived minimum ini- 
tial masses of 10 18 g for comet 109P/Swift-Tuttle and of 10 1B 
g for 55P/Tempel-Tuttle, based on the mass estimates for 
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the associated Perseids and Leonids meteor streams respec- 
tively, and assuming that the masses of the meteor streams 
were only a fraction of the masses still remaining in the re- 
spective nuclei. Asphaug and Benz (1996) modeled the tidal 
disruption of comet Shoemaker-Levy 9, and the formation 
of a long chain of clumps of fragments, during its near-miss 
with Jupiter in 1992. The authors assumed a stengthless 
body. By comparing the model results with observations, 
they obtained as a best-fit solution a diameter ~ 1.5 km 
and a density of ~ 0.6 g cm -3 for SW9. 

All the results mentioned before can only be applied 
to particular cases and are very uncertain. Since we do not 
know about binary comets or comet satellites, and no space- 
craft has measured the gravitational field of a comet (ex- 
cept for the Deep Impact experiment with comet 9P/Tempel 
1, as commented in section 6), the only general effect that 
can give us a direct estimate of the mass of the comet nu- 
cleus is the non-gravitational acceleration as shall be ana- 
lyzed in the next section. A method for computing cometary 
masses based on the non-gravitational effect was introduced 
by Rickman (1986, 1989), Rickman et al. (1987), and a 
somewhat similar approach was presented by Sagdeev et al. 
(1987). 

In this paper we plan to reevaluate the suitability of 
the non-gravitational force as a tool to estimate cometary 
masses following a procedure somewhat similar to that de- 
veloped by Rickman and colleagues. We will also analyze 
how the uncertainties involved in the different physical pa- 
rameters related to the problem affect the final results. Fi- 
nally, we will take advantage of the progress made in the 
determination of cometary sizes during the last decade to 
compute bulk densities from the derived non-gravitational 
masses. 

Section 2 presents a brief description of the non- 
gravitational forces. The method used, the model param- 
eters, and the data of the selected comets for our study are 
summarized in Section 3. In Section 4 we discuss the prob- 
lems in the cometary photomety and the empirical correla- 
tion between visual magnitudes and gas production rates. In 
Section 5 we present and discuss the visual lightcurves ob- 
tained for our sample of comets. The results are presented 
and discussed in Section 6, with a further analysis in Section 
7. Finally, a summary with our main conclusions is made in 
Section 8. 



2 NON-GRAVITATIONAL FORCES 

The main non-gravitational effect that can be detected in a 
periodic comet is a change in its orbital period, with respect 
to that derived from purely gravitational theory. Whipple 
(1950) proposed an icy conglomerate model for the cometary 
nuclei, and showed that the momentum transferred to the 
nucleus by the outgassing could cause the observed non- 
gravitational effect, an idea that was already advanced more 
than a century before by Bessel (1836). Because of the ther- 
mal inertia and the nucleus's rotation, Whipple showed that 
the direction of the net outgassing will be deviated a certain 
angle (known as the lag angle) with respect to the subsolar 
point. Today we know that some geometrical factors like the 
shape, the topography, the spin orientation or the location of 
the active areas on the nucleus surface, are also crucial in the 



determination of the lag angle, as suggested by thermophys- 
ical models (Davidsson and Gutierrez 2004, 2005, 2006, and 
Davidsson, Gutierrez and Rickman 2007). We emphasized 
that in this work the term 'non-gravitational forces' means 
the momentum transferred to the nucleus by the outgassing, 
since other non-gravitational forces acting on minor bodies 
of the Solar System are negligible in the case of comets. 

The change AP in the orbital period P due to the non- 
gravitational force can be expressed in terms of the radial J r 
and transverse Jt components of the non-gravitational accel- 
eration J by means of the Lagrangian planetary equations 
under the Gaussian form, which leads to 



AP 



J r sin(/)dt+ 



Jt 



dt 



,(1) 



where n is the mean motion, a the semimajor axis, e the 
eccentricity, / the true anomaly, and r the heliocentric dis- 
tance. The direction of the non-gravitational acceleration 
(which is in the opposite direction to that of the net out- 
gassing) can be described by the angle r\ with respect to the 
antisolar direction (i.e. the lag angle), and an azimuthal an- 
gle cj> in the plane perpendicular to the antisolar direction. 
The components of the non-gravitational acceleration can 
be expressed in terms of the angles 77 and <j> as follows 



J r — Jcos(rj), 

Jt — J sin(rj) cos((f>), 

Jn — J sin(rj) sin(0) 



(2) 



Equation (JTJ shows that AP depends on the radial and 
transverse components J r and Jt, but not on the normal 
component J n of the non-gravitational acceleration (i.e. the 
component perpendicular to the orbital plane). 

We find some confusion in the literature as regards to 
what we understand by AP. Some authors refer to a delay or 
advance in the time of the perihelion passage, instead of the 
change in the orbital period. These quantities are not equiv- 
alent since the change in the time of the perihelion passage 
is determined not only by the change in the period AP, but 
also by the change in the angular orbital parameters. 

Finally, according to Whipple (1950), we can relate the 
non-gravitational acceleration J and the comet's mass M by 
means of the conservation of momentum i.e. 

MJ = -Qmu, (3) 

where Q is the gas production rate, u is the effective outflow 
velocity, and m is the average molecular mass. 



3 THE METHOD 

In this section we introduce the assumptions we made in 
order to simplify the problem of the mass determination. 

By substituting J r and J t of eq. (JT| by the expressions 
given by eqs. @ and © , and then solving for the mass, we 
obtain the following expression 



M 



Qn\/1 — e 2 



a(l-e 2 ) J Q AP 



■ cos(?7) sin(/)di + 



Qmu 1 . , . , . 
-Tp-Bin(»/) cos(0)dt 



(4) 
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In order to solve the integrals in eq. (|4} we should know 
the outflow velocity and the angles as functions of time or 
the heliocentric distance. Our first approach will be to con- 
sider an average value < u > for all escaping molecules dur- 
ing an orbital revolution of period P. We also consider aver- 
age values for the angles < cos(?7) > and < sin(r;) cos(0) >. 
We assume these values as constants for all the comets of our 
sample. We also assume that water, as the principal volatile 
component in the cometary composition, dominates the gas 
production in the inner solar system (say r H 3 AU), so 
we take m as the water molecular mass. Therefore, we can 
express eq. Q as follows 

6n\/l — e 2 , rr , ,_. 

M ~ 2 A P m <U>(Ir + h), 5 

n AP 

where I r and It are given by 

e f p 
I r = , 2 , < cos(7;) > / Qsm(f)dt, 
a (t — e ) J 

f P Q 

It = < sin(ry) cos(0) > / — dt 
Jo r 

The computation of the effective outflow velocity < u > 
relies on sophisticated thermophysical modeling of the ex- 
change of linear momentum between the sublimating gas 
and the nucleus surface. The value of < u > depends not 
only on nucleus temperature, but also on shape and loca- 
tion of active areas, therefore is highly uncertain. As an ed- 
ucated guess, we consider < u > = (v t h, where vth is the 
mean thermal velocity (in an equilibrium gas), and £ is a 
dimensionless quantity that contains any difference between 
the physically meaningful parameter < u > and the arbi- 
trary chosen normalization quantity v t h- The mean thermal 
velocity is given by v t h = (8fcT/7rm) 1//2 , where k is Boltz- 
mann's constant, T is the temperature of the sublimating 
gases (which is about 200 K for the range of heliocentric 
distances of interest ~ 0.5 - 3 AU), and m, as before, the 
water molecular mass. For £ ~ 0.5 - 0.6 we obtain the value 
< u > ~ 0.25 km s _1 (Wallis and Macpherson 1981, Rick- 
man 1989). Some other works point to a somewhat higher 
value of < it > (around 0.3 km/s) (Peale 1989, Davidsson 
and Gutierrez 2004). Taking into account the range of pub- 
lished values, we finally adopted the mean value of (0.27 ± 
0.1) km s" 1 for the effective outflow velocity. 

As we said before, in general r\ and <f> will vary with 
the orbital position. Unfortunately, we cannot know the val- 
ues of these angles, unless we can study the comet nucleus in 
situ. Therefore, we depend on educated guesses and averages 
to sort this problem out. For most reasonable combinations 
of tentative values for the rotational period, the thermal in- 
ertia of the outermost layers of the nucleus and the heliocen- 
tric distance, Rickman (1986) and Rickman el al. (1987) ob- 
tained, by means of thermal models, rj ^ 30°. In more recent 
works with thermophysical models Davidsson and Gutierrez 
(2004, 2005, 2006) also found small values for t). These re- 
sults confirm our guess that the net non-gravitational force 
may be deviated from the radial direction just by a small 
angle. 

So we have only positive values for the parameter < 
cos(r;) >, while the parameter < sin(?7) cos(0) > has not 
such a constraint. This leads to an indetermination in the 
sign of the transverse non-gravitational contribution, as de- 



fined in eq. |(3J. To sort this problem out, we will use the 
following approach 

< sin(r;) cos(0) >= sig(AP) < | sin(^) cos(0)| > . 

To explain this assumption, we refer to eq. ([S}. There we 
can see that, for a comet with a symmetric curve Q(t) with 
respect to perihelion (this implies a symmetric lightcurve, as 
assumed in this work) the first integral will vanish, since a 
symmetric Q will also imply a symmetric J r , so the change 
AP will depend only on the transverse component J t . In 
this case it is easy to see that AP and Jt must have the 
same sign, and also AP and < sin(?7) cos(</>) >. 

Yet, as pointed out by Rickman (1986), most comet 
lightcurves are moderately or highly asymmetric (see for in- 
stance the lightcurves in Section 5), so the integral of the 
term containing J r will not vanish, and in some cases may 
become dominant. When one of the terms defined in eq. (JS| 
clearly dominates the other, there is no ambiguity in the re- 
sult: (i) when I r >> I t , the uncertainty in the sign of the 
transverse term contributes to the overall uncertainty in the 
mass as any other source of error, and (ii) when « I t , 
there is no ambiguity as explained before. The problem of 
the indetermination in the sign of the transverse term arises 
when both terms have similar absolute values, in which case 
we cannot easily discern if I t adds or subtracts to I r . An- 
other problem related to the case where I r « It arises be- 
cause the errors may make the difference I r - It either posi- 
tive or negative, leading for one of the signs to an unrealistic 
value M < 0. 

If we further assume a uniform random distribution in 
the range [0, 2-7r] for the angle cj>, we finally get the nominal 
values for our model parameters 

< cos(77) >^ 1, < j sin(77) cos(</>)| 0.1 

Now we turn our attention to the remaining parame- 
ters or quantities involved in the computation of the nu- 
cleus mass: the gas production curve, the relevant orbital 
parameters (semimajor axis and eccentricity), and the non- 
gravitational effect AP. Unlike the parameters discussed be- 
fore, these parameters are specific for each comet. 

As we can see from eq. ((5j), the computation of the 
cometary mass requires to know from observations the shape 
of the curve Q(t). There are only a few cases for which we 
have good measurements of gas production rates at the dif- 
ferent orbital positions, so in general we have to rely on 
lightcurves, and hence, assume some formula to approxi- 
mately convert total visual heliocentric magnitudes ran (i-e. 
the apparent magnitudes m corrected for the geocentric dis- 
tance A: mj, = m — 5 log A) to water production rates. In 
this regard we use the following empirical law, introduced 
by Festou (1986) 

log Q =aiXrrih + a-2 (6) 

We adopt for the coefficients a\ and 02 the values found by 
Jorda et al. (2008) (see Section 4.2 for a discussion of the 
empirical correlation between total visual magnitudes and 
gas production rates). In Table 1 we summarize the values 
adopted for the model parameters m, < u >, < \ cos(?7)| >, 

< I sin(ry) cos(</>)j >, a± and a^. 

The orbital elements of the selected comets are shown 
in Table 2. For the orbital parameters a and e we take the 
values corresponding to the epoch closer to the perihelion 
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passage. The symbols represent: r the time of the perihelion 
passage, q the perihelion distance (AU) , e the eccentricity, P 
the orbital period (years), uj the perihelion argument (°), Q 
the longitude of the ascending node (°), and i the inclination 
(°). The orbital elements and the epoch are referred to the 
2000.0 equinox. The elements were obtained from the cata- 
logue of cometary orbits of Marsden and Williams (2008). 
For those comets for which the lightcurves were based on 
photometric data corresponding to more than one appari- 
tion, an average of the values for each epoch was made. The 
mean motion n and the orbital period P were computed 
from the value used for the parameter a. 

The non-gravitational effect AP was computed from the 
non-gravitational orbital solution for the comet. In this case 
the equation of motion has three non-gravitational parame- 
ters Ai, A2, Az for the radial, transverse and normal compo- 
nents of the non-gravitational acceleration. The change AP 
will be independent of A3 and, for the symmetric case, also 
of Ai, so it will only depend on A2 through the equation 
(Marsden et al. 1973, Festou et al. 1990): 

AP = eWTEZ A r p a jr) dt 

n 2 J (j r 

The computed values of A2 are shown in Table 3. Al- 
though the model used by Marsden et al. (1973) (known as 
the standard symmetric model) was questioned in a physical 
sense to be not too realistic (since it assumes a symmetric 
outgassing curve with respect to perihelion), it has achieved 
in practice good fits to the astrometric observations of a 
large number of comets. 

The parameter A2 represents the transverse component of 
the non-gravitational acceleration at 1 AU from the Sun, 
and g(r) is an empirical function which describes the varia- 
tion of the water snow sublimation rate with respect to the 
heliocentric distance 
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(8) 



where a = 0.1113, m = 2.15, n = 5.093, k = 4.6142, and r 
= 2.808 AU, (Marsden et al. 1973). 

It may be argued that there may be a contradiction 
between the derivation of AP from the symmetric model 
and the use of a more realistic asymmetric model. The ex- 
planation of this seemingly contradictory procedure is that 
Ai, A2 arise from the best-fit orbital solution from Mars- 
den and Williams (2008) which is based on the symmet- 
ric model developed by Marsden et al. (1973). Therefore, 
for consistency AP has to be obtained from the A2 value 
for the symmetric model. We expect that AP values de- 
rived from asymmetric models should lead to a reasonable 
agreement with the values derived from eq. as far as 
we get good-quality orbit solutions for the comet's astro- 
metric positions based either on a symmetric or an asym- 
metric model. In order to check this assumption we also de- 
rived AP from Yeomans and Chodas's (1989) asymmetric 
model. In this regard the JPL Small-Body Database Browser 
(http://ssd.jpl.nasa.gov/sbdb.cgi) provides sets (Ai, Ai, 
Ai, DT) for a few Jupiter family comets derived from the 
latter model. DT represents the the time offset of the non- 
gravitational force maximum from the time of the perihelion 
passage. From the JPL web site we chose comets 6P, 10P, 
46P, 67P and 81P. The results (shown in Table 4) indicate 



that the values derived from both models are consistent, as 
expected. 



4 COMETARY PHOTOMETRY 

In this section we refer to the photometric data base and the 
criteria applied to construct the cometary lightcurves m/,(t). 
In this regard we analyze some problems in the determina- 
tion of cometary magnitudes, and discuss how we can relate 
these magnitudes to the gas production curve. 

4.1 Determination of cometary total magnitudes 

We understand as 'total' magnitudes, the ones that com- 
prise the light coming from the nuclear region and the 
extended coma. By contrast, the term 'nuclear' magnitude 
refers to the ones comprising only the central condensation. 
It should be noted that in most cases the latter magnitudes 
do not correspond to the bare nucleus, since it is very dif- 
ficult to avoid contamination by coma light. Nevertheless, 
efforts are being made to determine true magnitudes of 
comet nuclei by observing them far from the Sun or by a 
coma substraction method (e.g. Tancredi et al. 2006). 

The precise determination of cometary total magni- 
tudes presents serious problems since comets do not ap- 
pear as point sources like the stars but as nebulous sources. 
The measured magnitude will depend on the instrument em- 
ployed (e.g. Bobrovnikov 1941, Morris 1973). A telescope of 
small aperture will bring within the field of view the central 
condensation of the comet plus the extended coma, while a 
telescope of large aperture will only bring the central con- 
densation. In the latter case the comet will thus look fainter 
than in the former case. This aperture effect varies both from 
comet to comet and among different telescopes types. The 
average aperture effect also depends upon the degree of con- 
densation of the comet. Therefore magnitudes measured by 
different observers with different instruments can differ by 
several units of magnitude. As an example, Fig. 1 shows the 
apparent total magnitudes for the 2005 apparition of comet 
9P/Tempel 1, measured by a visual method (i.e. an esti- 
mate made with a small telescope or binoculars by eye) as 
a function of the instrumental aperture. All the magnitude 
data used in this work comes from the International Comet 
Quaterly (ICQ) (courtesy of D. Green), unless we specify 
another source. Following the ICQ's recommendations, to 
minimize the aperture effects the magnitudes should be ob- 
tained using the smallest objective aperture and magnifica- 
tion needed to easily see the comet. In Fig. 1 we can see a 
general trend to reporting brighter magnitudes with smaller 
apertures, and viceversa. We also note that even for a given 
aperture, there is a noticeable dispersion in the magnitudes 
reported by different observers. 

We are restricted to the spectral region where the hu- 
man eye is more sensitive, i.e. all the magnitudes considered 
in this work are visual magnitudes. So in the following, by 
'visual magnitudes' we refer to the cometary magnitudes es- 
timated by a visual method, while by 'CCD magnitudes' we 
refer to estimates based on a CCD detector with a V broad- 
band filter or on an unfiltered CCD detector. According to 
the ICQ, it has been found by some observers that V and 



Cometary masses derived from non- gravitational forces 5 



Table 1. Summary of model parameters. 



Description Symbol Nominal value Unit 

Water molecular mass m 2.98897 x 1CT 26 kg 

Effective outflow velocity (averaged value) < u > 0.27 km s _1 

Angular parameter < cos(?7) > 1.0 

Angular parameter < | sin(rj) cos(</>)\ > 0.1 

Calibration coefficient ai -0.2453 

Calibration coefficient 30.675 



Table 2. Orbital parameters of the studied comets. 



Comet 


T 


(I 


e 


P 




n 


i 


Epoch 


ip 


1986 Feb. 9.4589 


0.587104 


0.967277 


76.0 


111.8657 


58.8601 


162.2422 


86 Feb. 19 


2P 


1990 Oct. 28.5678 
1994 Feb. 9.4785 
1997 May 23.5979 
2000 Sept. 9.6669 


0.330890 
0.330918 
0.331397 
0.339539 


0.850220 
0.850211 
0.850013 
0.846898 


3.28 
3.28 
3.28 
3.30 


186.2348 
186.2723 
186.2733 
186.4842 


334.7492 
334.7273 
334.7202 
334.5993 


11.9450 
11.9401 
11.9293 
11.7554 


90 Nov. 5 
94 Feb. 17 
97 June 1 
00 Sept. 13 


6P 


1995 July 27.3234 


1.345813 


0.614042 


6.51 


178.0502 


138.9895 


19.5239 


95 July 22 


9P 


1994 July 3.3141 
2005 July 5.3145 


1.494151 
1.506166 


0.520255 
0.517491 


5.50 
5.52 


178.9009 
178.8382 


68.9864 
68.9380 


10.5519 
10.5301 


94 June 17 
05 July 9 


10P 


1983 June 1.5367 
1988 Sept. 16.7345 
1999 Sept. 8.4214 


1.381401 
1.383426 
1.481678 


0.544895 
0.544429 
0.522818 


5.29 
5.29 
5.47 


190.9473 
191.0643 
195.0232 


119.8323 
119.7925 
118.2114 


12.4338 
12.4282 
11.9767 


83 May 26 
88 Oct. 6 
99 Sept. 19 


19P 


1987 Dec. 18.3239 
1994 Nov. 1.4946 
2001 Sept. 14.7306 


1.356795 
1.365126 
1.358204 


0.624213 
0.622800 
0.623891 


6.86 
6.88 
6.86 


353.3370 
353.3587 
353.3750 


75.4332 
75.4238 
75.4249 


30.3254 
30.2708 
30.3248 


87 Dec. 31 
94 Oct. 15 
01 Sept. 8 


22P 


1983 Aug. 10.3011 
1996 July 2.1914 


1.576327 
1.579573 


0.544526 
0.544069 


6.44 
6.45 


162.8877 
162.8377 


120.9308 
120.9065 


4.7210 
4.7211 


83 Aug. 14 
96 July 16 


45P 


1990 Sept. 12.6843 
1995 Dec. 25.9844 
2001 Mar. 29.9269 


0.541252 
0.532050 
0.528413 


0.821869 
0.824231 
0.825080 


5.30 
5.27 
5.25 


325.7863 
326.0532 
326.1311 


89.3080 
89.1561 
89.0828 


4.2198 
4.2484 
4.2556 


90 Sept.26 
95 Dec. 29 
01 Apr. 1 


46P 


1991 Sept. 20.6255 
1997 Mar. 14.1499 


1.083306 
1.063769 


0.652242 
0.656748 


5.50 
5.46 


356.1687 
356.3418 


82.2930 
82.2051 


11.6820 
11.7226 


91 Sept.21 
97 Mar. 13 


67P 


1982 Nov. 12.0993 
1996 Jan. 17.6560 


1.306144 
1.300032 


0.629152 
0.630193 


6.61 
6.59 


11.3667 
11.3861 


51.0145 
51.0070 


7.1170 
7.1135 


82 Nov. 7 
95 Dec. 29 


81P 


1990 Dec. 16.9164 
1997 May 6.6277 
2003 Sept.25.9302 


1.578061 
1.582623 
1.590367 


0.540980 
0.540220 
0.538787 


6.37 
6.39 
6.40 


41.6421 
41.7681 
41.7516 


136.2048 
136.1564 
136.1413 


3.2437 
3.2424 
3.2400 


90 Dec. 15 
97 Apr. 22 
03 Oct. 8 



unfiltered magnitudes of comets do not differ by more than 
several tenths of a magnitude. Mikuz and Dintinjana (2001) 
have compared also CCD V and CCD unfiltered observa- 
tions finding no meaningful differences among them. 

For faint comets, detectors (e.g. CCDs) tend to record 
only the nuclear condensation while losing the broad (faint) 
coma that fades into the sky background. Therefore, the 
measured 'total' magnitudes will be fainter than they actu- 
ally are. The question of how CCD magnitudes - whether fil- 
tered or unfiltered - relate to visual magnitudes is still open. 
Towards the limit of visual observations (near magnitudes 12 
- 15), it has been noted that CCD total magnitudes are typ- 
ically 1-3 magnitudes fainter than visual estimates (Green 
1996). As an example, Fig. 2a shows the apparent visual 
magnitudes as compared to the apparent CCD magnitudes, 



for the 2005 apparition of 9P/Tempel 1. The figure shows 
that the CCD magnitudes, as a group, are systematically 
fainter than the visual magnitudes. The CCD magnitudes 
are also more scattered, while the visual magnitudes depict 
a more compact group. 

To take into account the sky fading effect, Kresak and 
Kresakova (1989) suggest the following empirical correction 
for the measured apparent magnitudes of faint comets with 
m > 9 

0.5m + 4.5 (9) 

Figure 2b shows that the corrected CCD magnitudes are 
much less dispersed (resulting in a more flat curve). Yet, in 
spite of a good agreement with respect to the visual magni- 
tudes group, as we go backwards in time ( > 50 days before 
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Table 3. Non-gravitational effect in the studied comets. 





A T\r\ ti + i r\Y\ 
XV L/dl 1 LHJI1 


A 2 xio 8 a 


AP 






(All/days 2 ) 


(days) 


lP/Halley 


1986 


+0.0155 


+4.1063 


2P/Enckc 


1990 


-0.0016 


-0.0048 




1994 


-0.0007 


-0.0021 




1997 


-0.0007 


-0.0021 




9000 


n 0007 


n nri9i 

-U.UUZ 1 


6P/d' Arrest 


1995 


+0.0993 


+0.1297 


9P/Tempel 1 


1994 


+0.0017 


+0.0013 




2005 


+0.0017 


+0.0013 


lOP/Tempel 2 


1983 


+0.0014 


+0.0013 




1988 


+0.0014 


+0.0013 




1999 


+0.0014 


+0.0011 


19P/Borrelly 


1987 


-0.0376 


-0.0522 




1994 


-0.0473 


-0.0651 




2001 


-0.0473 


Ofi^ 


22P/Kopff 


1983 


-0.1127 


-0.0962 




1996 


-0.1073 


-0.0912 


45P/H-M-P 


1990 


-0.0553 


-0.2042 




1995 


-0.0505 


-0.1887 




2001 


-0.0505 


-0.1893 


46P/Wirtanen 


1991 


-0.1673 


-0.2557 




1997 


-0.1373 


-0.2137 


67P/C-G 


1982 


+0.0096 


+0.0137 




1996 


+0.0096 


+0.0138 


81P/Wild 2 


1990 


+0.0092 


+0.0077 




1997 


+0.0092 


+0.0077 




2003 


+0.0092 


+0.0076 



a Source: Marsden and Williams (2008). 
Table 4. Non-gravitational effect according to two different models. 



Comet Symmetrical model Asymmetrical model 







(Marsden et al 1973) 






(Yeomans & Chodas 1989) 






Apparit: 


ion Parameters' 1 


AP 


Apparition 




Parameters 6 




AP 






A 2 xlO 8 






Ai x 10* 


A 2 x 10 8 


DT 








(AU/days 2 ) 


(days) 




(AU/days 2 ) 


(days) 


(days) 


6P/d' Arrest 


2008 


+0.0993 


+0.1288 


2008 


+0.3804 


-0.0117 


+97.3 


+0.1302 


lOP/Tempel 2 


2005 


+0.0014 


+0.0012 


2006 


+0.0325 


-0.0010 


+24.3 


+0.0011 


46P/Wirtanen 


2008 


-0.1373 


-0.2148 


2007 


+0.3735 


-0.1302 


-4.4 


-0.2159 


67P/C-G 


2009 


+0.0096 


+0.0145 


2008 


+0.1160 


-0.0062 


+38.8 


+0.0171 


81P/Wild 2 


2003 


+0.0092 


+0.0076 


2006 


+0.1489 


+0.0333 


-60.6 


+0.0078 



Source: Marsden and Williams (2008). 
Source: JPL Small-Body Database Browser. 
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perihelion) , we see an increasing departure of the CCD group 
from the visual group: the CCD lightcurve pre-perihelion 
branch seems to be significantly different from the visual 
pre-perihelion branch, with a less steep slope. We conclude 
that the corrected CCD lightcurve does not agree well with 
respect to the visual lightcurve. 

Since we could not find an adequate correction to CCD 
magnitudes for the sky fading effect, we use only the vi- 
sual observations in this work. We also found that for the 
comets of our sample the visual observations are much more 
numerous than the CCD observations, during the most ac- 
tive period of the comet (see Fig. 2a for instance). 

But the visual observations also show a considerable 
dispersion, as we saw before (the observations are made 
with different instrumentation, methodology and/or refer- 
ence stars, and probably in different local conditions) , so we 
have to adopt some criteria to homogenize the visual mag- 
nitude data. 

As the first step, we discarded the observations made 
when the comet was at low altitude without an atmospheric 
extinction correction applied, or when the magnitude esti- 
mate is not very accurate, or when it was made under poor 
weather conditions according to the information provided 
by the observer. 

Whereas there are not obvious effects that may lead to 
an increase in the perceived brightness, other than observa- 
tion during an outburst, there are several effects (like large 
apertures or an excess of magnification as we mentioned be- 
fore, and others like moonlight, twilight, haze, cirrus clouds, 
dirty optics, etc.) that will decrease the brightness and hence 
overestimate the observed visual magnitudes. As noted by 
Ferrm (2005), the top observations follow a smooth trend, 
defining a sharp boundary, while the lower part of the visual 
observations is more diffuse, poorly defined (see the visual 
lightcurve of Tempel f in Fig. 2a). Following Rickman et 
al. (1987), we take a third, fourth or sixth order polynomial 
fit to the upper envelope of the ensemble of photometric 
observations as the lightcurve m^(t). In practice this poly- 
nomial is obtained as a least-square fit to a set of selected 
data points. This set is defined by the following procedure: 
we divide the time domain (defined by the observations) in 
N bins of a fixed size. Then we select, for each one of these 
time bins, the three brightest observations (namely the ones 
with the lowest values of the heliocentric magnitude), that 
do not depart more than 3-sigma from the mean value of the 
observed magnitudes within the bin. We neglect those time 
bins with very few observations. Therefore, we allow for a 
certain range of magnitudes to define our upper envelope, 
not just the brightest one within each time bin, which gives 
us more confidence that such a range should cover all po- 
tential sources of errors in the magnitude measurements. We 
used a time bin of 3 days, except for comet 46P/Wirtanen 
for which we used a time bin of 5 days due to the relatively 
small number of observations involved. 

The definition of the lightcurve, as a third, fourth or 
sixth order polynomial, implies that it varies smoothly with 
time, so comets suspected of having outbursts or another 
" anomalous" photometric behavior are excluded of our sam- 
ple. The polynomial fit does not have any symmetry restric- 
tion either, so asymmetric fits (with respect to perihelion) 
are allowed, and actually are usually the case. The only in- 
convenience we find with this lightcurve definition is that we 



can not extrapolate the magnitude beyond the time interval 
defined by the observations, which is equivalent to neglect 
the gas production outside such an interval. Nevertheless, 
outside the observational time interval the gas production 
usually drops several orders of magnitude with respect to 
the maximum, so we do not expect to introduce a significant 
bias in the computed mass due to this constraint (in Sec- 
tion 6 we discuss the uncertainty in the computed cometary 
masses due to the assumptions we make as regards to their 
lightcurves) . 

We also neglect short-term variations in the visual mag- 
nitudes (like the amplitude of the rotational lightcurve), and 
variations due to the phase angle, since these sources of in- 
trinsic brightness scattering would be included in the fitted 
lightcurve by considering sets of the three brightest magni- 
tudes in every time bin, as discussed before. 

4.2 Visual magnitudes and gas production rates 

As we pointed out before, at heliocentric distances r H 3 
AU the cometary activity is governed by the sublimation 
of water ice. The gas production rate of water ice is thus 
an adequate quantitative indicator of cometary activity. Di- 
rect observations of the H2O production rate are difficult 
and sparse because they involved sophisticated and over- 
subscribed instrumentation. Indirect observations (e.g. ob- 
servation of water-derived products like the OH radical from 
its radio lines at 18 cm and from narrow-band photometry 
in the near UV) are more accessible, but also limited. The 
huge database of visual magnitudes gives another approach 
to this problem (Crovisier 2005). In this regard many au- 
thors have studied the empirical relation between total vi- 
sual magnitudes (corrected for geocentric distance) and gas 
production rates (Festou 1986, Roettger et al. 1990, Jorda 
el al. 1992 and 2008, de Almeida et al. 1997). From a sta- 
tistical analysis of 37 comets, based upon values from 
the ICQ and Qh 2 o from the Nangay database, Jorda et al. 
(2008) derived the linear relation shown in eq. ((6]l where 
Q — Qh 2 o, a-i = -0.2453 and 02 = 30.675. This correla- 
tion is based on a larger database than that of the former 
work of Jorda et al. (1992), and it would be more reliable 
(L. Jorda, private communication). This relation (or simi- 
lar ones) has been used to predict gas production rates of 
comets in the absence of direct measurements (e.g. Festou et 
al. 1990, Jorda and Rickman 1995, de Almeida et al. 1997). 
It has also been explored a similar relationship for the CO 
molecule by Biver (2001), for r > 3 AU. 

Festou (1986) argued that the relationship (eq.|SJ seems 
to hold very well if a set of conditions is fulfilled: 1) a steady 
state is established in the coma (which would imply that the 
lightcurve varies smoothly with time), 2) the visible part of 
the comet spectrum is dominated by C2 emissions, and 3) 
the ratio of C2 to OH does not vary from comet to comet and 
with the heliocentric distance. This latter condition could 
suggest a similar composition for all comets. In this regard 
A'Hearn et al. (1995) found that there is little variation of 
relative abundances with heliocentric distance, and there is 
also little variation from one apparition to the next for most 
short-period comets. 

In Fig. 3 we show the empirical water production rate 
Q(t) as a function of time that we obtained for comets 
9P/Tempel 1, 19P/Borrelly, 67P/C-G and 81P/Wild 2. 
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Figure 1. Apparent visual total magnitudes (circles) as a function of time from the perihelion passage for 
the 2005 apparition of 9P/Tempel 1 (data from the ICQ archive). The size of the circles indicates the size of 
the aperture used: the largest circles corresponds to apertures of 40 $J diameter < 45 cm, while the smallest 
circles corresponds to apertures of diameter < 10 cm. 
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Figure 2. Apparent visual total (asterisks) and CCD (plus signs) magnitudes as a function of time (days 
from the perihelion passage), for the 2005 apparition of 9P/Tempel 1 (data from the ICQ archive), (b): 
The same as the previous figure, except that the CCD magnitudes were corrected for the sky fading effect 
according to Kresak and Kresakova (1989)'s formula. 
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Each curve Q(t) was obtained from the respective visual 
lightcurve (shown in Section 5), where the heliocentric vi- 
sual magnitudes (rrih) were converted to water production 
rates through eq. JBJ, using Jorda et al. (2008) calibration's 
coefficients. Fig. 3 also shows an individual calibration for 
each comet, obtained from a least-square linear fit between 
the observed log(Q_ff 2 o) and the estimated magnitudes rrih 
at the respective observing times (we did not applied any fil- 
ter or weight to the compiled observational data). We cannot 
tell at this point if these individual calibrations make a bet- 
ter fit than the Jorda et al. (2008) calibration, due to the 
scatter (and relatively scarcity) of the observational mea- 
surements. We also computed the water production rates 
predicted by the Jorda et al. (1992) calibration's coefficients, 
and found that the residuals between predicted and observed 
water production rates were about a factor of two larger than 
the residuals obtained with the Jorda et al. (2008) constants. 
So we chose the Jorda et al. (2008) calibration as our nom- 
inal one, for all comets of our sample, and take the differ- 
ences in the computed masses (with respect to the computed 
masses using the individual calibrations), as an estimate of 
the uncertainty due to the calibration's coefficients (see Sec- 
tion 6.1 for more details on the error's estimate). 

We consider, as Crovisier (2005), that the empirical re- 
lation (eq. [6]), even if it is physically poorly understood, it 
is a very useful tool. Maybe the present lack of a physi- 
cal interpretation of this empirical relationship reflects our 
poor knowledge of how the cometary activity takes place. 
By assuming a very simplified physical model for the comet 
nucleus, that consists of a spherical nucleus with a smooth 
surface, and that most of the solar radiation is evenly spread 
across the surface and spent in the sublimation of ices, we 
then have Q oc r~ 2 (r H 3 AU), so logQ oc -21ogr. Since 
m,h — H + 2.5nlog r, where H is the absolute total magni- 
tude and n is the photometric index that gives the comet's 
brightness variation with r, we can relate the previous two 
relations obtaining 

Q g 

logQ xm h + C, (10) 

n 

where C is independent of r and only depends on H. Eq. 
(| 10p allows us to relate the index n with the coefficient ai of 
the empirical calibration (eq. |BJ|. It is interesting to stress a 
few consequences of this very simple approach: if the comet 
brightens exactly like a bare solid body (i.e. the solid nucleus 
without gaseous activity), we have n = 2, and then a± — -0.4. 
But comets usually show activity indices n > 2 which indi- 
cates that when approaching the Sun they brighten much 
more than expected for a bare solid body, because in this 
case most of the comet light comes from sunlight scattered 
by dust particles in the coma and fluorescent emission from 
gaseous molecules, which greatly increases their photometric 
cross-section. So we expect n > 2 and then ai > -0.4. This 
is consistent with the value ai — -0.2453 from Jorda et al. 
(2008), for which we derive n ~ 3. This latter value is closer 
to the average of the observed slopes of a large number of 
comets lightcurves (n — 4). As it was pointed out by Biver 
(2001), the ultimate goal would be to find the scientific jus- 
tification of this photometric index. This should be related 
to physical parameters like the amount of dust relative to 
gas, the size distribution of the dust particles, the relative 
abundance of C2, etc, as a function of the heliocentric dis- 



tance. Of course all these implies a much more sophisticated 
model for the comet nucleus that is beyond the scope of this 
work. 



5 THE STUDIED COMETS AND THEIR 
LIGHTCURVES 

In this section we present the criteria applied to select the 
comets for this study, and the lightcurves we obtained for 
the selected comets. 

We will limit this study to the short-period comets (i.e. 
P < 200 yr) since these comets have been observed in more 
than one apparition, so for them it is possible to obtain 
a good measurement of the non-gravitational effect APl]. 
To select the comets, the following set of conditions was 
adopted: 1) a measured non-gravitational change AP, 2) 
an adequate photometric coverage during the most active 
period around the perihelion passage, 3) a good estimate 
of nucleus size (in order to compute a bulk mass density), 
and preferentially (but not exclusive) 4) a measured water 
production rate (in order to check the coefficients for 
the empirical calibration of eq. |BJ. The comets selected 
were lP/Halley, 2P/Encke, 6P/d'Arrest, 9P/Tempel 1, 
lOP/Tempel 2, 19P/Borrelly, 22P/Kopff, 45P/Honda- 
Mkros-Pajdusakova, 46P/Wirtanen, 67P/Churyumov- 
Gerasimenko and 81P/Wild 2. 

The database of apparent total magnitudes used in this 
work is from the ICQ archive for observations after 1990, 
and from the Comet Light Curve Catalogue/ Atlas Part. I. 
(Kamel 1991a) for observations prior to 1990. 

For some comets we use observations involved in more 
than one apparition to make a composite lightcurve, to com- 
pensate for the relative low number of observations in a sin- 
gle apparition. This is possible since no significant variation 
of the orbital elements is noticed on two or more consec- 
utive observed apparitions of the chosen comets, and also 
by assuming that the shape of the lightcurve remains al- 
most unchanged during consecutive perihelion passages of 
the comet. 

Figures 4-14 show the lightcurves obtained as plots 
of the total visual heliocentric magnitude mj, as a function 
of the time relative to the perihelion passage t (the obser- 
vations were previously filtered according to the criteria ex- 
plained in Section 4). The horizontal line indicates the cut- 
off magnitude rrihc (i.e. the polynomial fit is considered a 
good approximation to the lightcurve only within the inter- 
val where rrih ^ rrihc)- For the conversion from apparent to 
heliocentric magnitudes we compute the comets' distances 
with the Mercury orbital integrator (Chambers 1999), and 
use some routines from Press et al. (1992) for function in- 
terpolation, and other routines from the NOVAS package of 
the U.S. Naval Observatory, written by G.H. Kaplan, for the 
conversion of time from Gregorian to Julian date. 

1 The evaluation of non-gravitational forces is more difficult for 
long-period comets, since these have not been observed in a sec- 
ond apparition to check for a change in the orbital period. Never- 
theless, non-gravitational terms have been fitted to the equations 
of motion of several LP comets leading to more satisfactory or- 
bital solutions. 
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Figure 3. Empirical water production rate as a function of time Q(t) (full line), superposed to observational 
data (small circles), for the 2005 apparition of 9P/Tempel 1 (top, left), for the 1987, 1994 and 2001 
apparitions of 19P/Borrelly (top, right), for the 1982 and 1996 apparitions of 67P/C-G (bottom, left), 
and for the 1990, 1997 and 2003 apparitions of 81P/Wild 2 (bottom, right). The empirical function Q(t) 
was obtained from the visual lightcurve (shown in Section 5) via eq. JB), using Jorda et al. (2008) calibration's 
coefficients. The dotted line shows an individual calibration for each comet, obtained from a least-square 
linear fit between the observed log(Q h 2 o) an d the estimated m/, from the lightcurve at the same observing 
times. The data corresponds to observed water production rates extracted from Davidsson et al. (2007) for 
comet Tempel 1 (courtesy of P. J. Gutierrez), from Davidsson and Gutierrez (2004) for comet Borrelly, from 
Davidsson and Gutierrez (2005) for comet C-G, and from Davidsson and Gutierrez (2006) for comet Wild 
2. 



As we discussed in Section 3, most comet lightcurves 
are asymmetric with respect to the perihelion. As a first ap- 
proach, the perihelion asymmetries of the comet lightcurves 
could be mainly related to structural changes of the nu- 
cleus that occur during the approach to the Sun, like the 
expulsion of an insulating dust mantle that would leave 
exposed areas of fresh ice. This effect, combined with the 
thermal inertia, would increase the gaseous activity after 
the perihelion passage. Most of the studied comets effec- 
tively show a post-perihelion brightness excess. These are 
the cases of comets lP/Halley, 6P/d Arrest, lOP/Tempel 
2, 19P/Borrelly, 45/H-M-P, and 67P/C-G. On the other 
hand, comets 2P/Encke, 22P/Kopff and 81P/Wild 2 show 
a pre-perihelion brightness excess. Hence the gaseous activ- 
ity could have its maximum before perihelion. This could be 
taken as another evidence of the complexity of the cometary 
activity, suggesting that other causes, like the concentration 
of the gaseous sublimation in a few discrete active areas, 
and the orientation of the spin axis, which may, or may not, 
favor the illumination of the active areas, may also play a 
role in the asymmetric behavior. Comets 9P/Tempel 1 and 
46P/Wirtanen show the most symmetric curve of the sam- 
ple. Comet Tempel 1 shows a slightly pre-perihelion excess 
of activity, but in this case the near- and post-perihelion 



coverage is very poor, due to unfavorable geometrical con- 
ditions. 

The maxima of the comet lightcurves seem to occur typ- 
ically around or about a few days before or after perihelion, 
although greater departures are common (see for example 
the lightcurves of comets d' Arrest and C-G), as it was al- 
ready noticed by Festou (1986). 

We note that comet Halley at its maximum brightness 
is several magnitudes brighter than the rest of the comets of 
the sample at their maxima. This should be correlated with 
a much greater gaseous activity for Halley as compared to 
the other comets. For instance, Halley shows a gas produc- 
tion rate of about 10 31 s -1 near perihelion, while the rest 
of the studied comets are approximately in the range [0.1 - 
2.0] xlO 28 s _1 near perihelion (A'Hearn et al. 1995, Crovisier 
et al. 2002). This result may be partly due to the different 
comet sizes, and partly due to the different source regions: 
while Halley is the prototype of the population of the H alley- 
Type (HT) comets - for which the source region would be the 
Oort Cloud (Fernandez 2005, Ch. 7), the rest of the studied 
comets belongs to the Jupiter Family (JF) population - for 
which the main source region would be the trans-neptunian 
belt (Fernandez 1980, Duncan et al. 1988). Since JF comets 
may have spent up to hundreds or thousands of revolutions 
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on short-period orbits near the Sun, they may be more phys- 
ically evolved than HT comets, so lP/Halley may still be a 
comet with a fresher and more active surface, leading to its 
observed substantially higher activity. 

Because of their short orbital periods (P < 20 yr), most 
of the discovered JF comets have been observed in several 
apparitions, and in some cases through all the orbit up to 
their aphelia, allowing to gather a very valuable wealth of 
physical data. Furthermore, several JF comets have been 
or will be the targets of space missions, which will help to 
greatly increase our knowledge about their physical nature. 
This is the case of four out of the eleven comets studied in 
this work. In the following we make a brief summary for 
each one of the lightcurves obtained. 



5.1 lP/Halley 

Comet lP/Halley was the target of several international 
space missions during its last apparition in 1986, which 
could take images of its nucleus of irregular shape, and con- 
firmed, in general terms, Whipple's model. From the ground 
the comet was bright enough to be visible to the naked 
eye, though far below the spectacular previous apparition 
of 1910. This comet presents the best photometric cover- 
age in one single apparition, as compared to the rest of the 
studied comets, resulting in a good quality lightcurve. The 
plot in Fig. 4 contains 1221 observations corresponding to 
the 1986 apparition, and shows the fitted lightcurve m,h{t). 
The lightcurve obtained is similar to others found in the 
literature (e.g. Ferrin 2005). 



5.2 2P/Encke 

For this comet we combined the observations of four con- 
secutive apparitions: 1990 (186 observations), 1994 (202 ob- 
servations), 1997 (29 observations), 2000 (24 observations). 
The observations of 1997, in spite of their low number, were 
important to define better the post-perihelion branch, for 
which only the 1994 apparition has post-perihelion obser- 
vations. Figure 5 shows a good consistency between the 
different apparitions, confirming our assumption that the 
lightcurve has not changed significantly from one apparition 
to the next, during the last decade. Yet we observe a secular 
decrease of the non-gravitational effect during this period. 
This was already noticed by Kamel (1991b) who studied 
the lightcurve evolution of this comet from 1832 to 1987. 
Kamel also found a secular decrease of the perihelion asym- 
metry (which he related to a shift in the time of maximum 
brightness from 3 weeks before perihelion to a few days after 
perihelion in the latest apparitions). He concluded that the 
shape of the lightcurve has not significantly changed during 
short intervals of time (like the last apparitions he studied), 
in particular the perihelion asymmetry has remained, show- 
ing a faster brightness decrease after the maximum and a 
slower increase before it. Our composite lightcurve, based on 
most recent observational data, is in good agreement with 
the results from Kamel. 



5.3 6P/d' Arrest 

The lightcurve dataset for this comet contains 424 obser- 
vations from the 1995 apparition (the last apparition in 
2002 was unfavorable since the comet was in conjunction 
with the Sun close to the time of the perihelion passage). 
The lightcurve is shown in Fig. 6. This comet shows an ex- 
traordinary strong asymmetry with respect to perihelion. 
The brightness rises rapidly as the comet passes perihelion, 
reaching the maximum around ~ 40 days after perihelion, 
with a slow decline after the maximum. The lightcurve ob- 
tained is similar to others found in the literature (e.g. Szu- 
towicz and Rickman 2006). 

5.4 9P/Tempel 1 

This comet has been the most recent target of a success- 
ful space mission, the Deep Impact, which encountered the 
comet about 1 day before perihelion on July 4, 2005. There- 
fore a great amount of data has been collected for this comet, 
including a good photometric coverage. Unfortunately, like 
the former apparition of 1994, there are very few post- 
perihelion observations due to geometrical circumstances. 
The lightcurve dataset for this comet contains 447 observa- 
tions from the 2005 apparition. The lightcurve is shown in 
Fig. 7. We also studied the lightcurve corresponding to the 
1994 apparition (with a dataset of 603 observations), and 
found an average brightness decrease of about 1-2 mag from 
1994 to 2005. The perihelion asymmetry seems to remain, 
with a slightly pre-perihelion predominance, but this feature 
is quite uncertain since there are not enough post-perihelion 
observations to define a good post-perihelion lightcurve. 
Ferrin (2005, 2007) has published a lightcurve for this comet, 
but only based on the 1994 data, which is similar to the 1994 
lightcurve we obtained. 

5.5 lOP/Tempel 2 

The lightcurve dataset for this comet contains 153 observa- 
tions from the 1983 apparition and 143 observations from the 
1988 apparition. The composite lightcurve is shown in Fig. 8. 
There are few pre-perihelion observations in 1983, but those 
from 1988 are enough to properly define the pre-perihelion 
branch of the lightcurve. We found a good consistency be- 
tween the observations of both apparitions, which implies 
that the lightcurve remains practically the same. Rickman 
et al. (1991a) have studied the lightcurve of this comet dur- 
ing 13 consecutive apparitions from 1899 up to 1988. They 
find that after perihelion the 1983 and 1988 apparitions yield 
quite different lightcurves. The disagreement between their 
result and ours may be due to a different filtering criteria 
applied to the observations (since the original database - 
Kamel 1991a - is the same). In correspondence with a neg- 
ligible variation of the lightcurve between both apparitions, 
we do not find any significant change in the orbital param- 
eters nor in the non-gravitational effect (see Tables 2 and 
3, respectively). Anyway, both works agree in the remark- 
able asymmetry showed by this comet lightcurve, which is 
very similar to the dArrest case: a rapid brightness increase 
before the maximum, followed by a slow decrease after it. 

We also studied the 1999 lightcurve, based on 239 ob- 
servations (the 1994 and 2005 apparitions were too poor in 
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number of observations). We found a significant brightness 
decrease of 1-2 mag from 1988 to 1999, and a shift in the 
time of the maximum brightness from around 14 days after 
perihelion to a few days after perihelion. Notwithstanding 
this remarkable change in the lightcurve, the asymmetry re- 
mains in shape: faster brightness increase before the max- 
imum, and a slower brightness decrease after it (although 
the pre-perihelion increase is less steep than that of the 1988 
lightcurve). This brightness decrease could be related to a 
small but noticeable increase (of about 7 %) of the perihelion 
distance from 1988 to 1999 (see Table 2). We also noticed 
a certain decrease (of about 20 %) of the non-gravitational 
change AP (see Table 3). This comet shows the smallest 
non-gravitational effect of the sample (closely followed by 
Tempel 1). Rickman et al. (1991a) suggested a possible cor- 
relation between the lower values of the non-gravitational 
effect on Tempel 2 and the apparent drop in the gas pro- 
duction as inferred from the lightcurves. 

We conclude that there is a remarkable difference of the 
1983 and 1988 composite lightcurve with respect to the 1999 
lightcurve. Since the quality of the former is much better, 
we adopt it as our nominal lightcurve (a short explanation 
of how we evaluate the lightcurve quality is given at the end 
of this section). 

5.6 19P/Borrelly 

This comet has been the target of the Deep Space 1 mis- 
sion, which encountered the comet on September 22, 2001, 
about 8 days after perihelion. The lightcurve dataset for this 
comet contains 360 observations from the 1987 apparition, 
853 observations from the 1994 apparition, and 242 obser- 
vations from the 2001 apparition. The composite lightcurve 
is shown in Fig. 9. We found no meaningful changes in the 
lightcurve from one apparition to the next, during the stud- 
ied period. We discarded a small group of bright measure- 
ments of the 2001 apparition bunched together around 20 
days before perihelion as a sort of spike, since they depart 
from the general trend of the bulk of the observations. This 
spike could be due to a little outburst produced around 30 
to 20 days before perihelion. Borrelly shows the same asym- 
metrical behavior as, for example, d'Arrest or Tempel 2: a 
faster brightness increase before the maximum followed by a 
slower decrease after it. The composite lightcurve obtained 
in this work resembles very much the one from Ferrin (2005) . 

5.7 22P/Kopff 

The lightcurve dataset for this comet contains 349 observa- 
tions from the 1983 apparition and 586 observations from 
the 1996 apparition. The composite lightcurve is shown in 
Fig. 10. This comet presents a curious lightcurve since, even 
though it shows a faster brightness increase before the max- 
imum, like those comets with a post-perihelion excess (e.g. 
C-G, d'Arrest, Tempel 2 or Borrelly), it shows however a 
pre-perihelion excess. We note however that the brightness 
maximum is quite uncertain, since the comet seems to hold 
a similar activity level for several weeks, so we cannot ex- 
clude the possibility that the real maximum could be a few 
days after perihelion. Anyway, it is clear from Fig. 10 that 
the comet shows more activity before than after perihelion. 



5.8 45P/Honda-Mkros-Pajdusakova 

The lightcurve dataset for this comet contains 117 observa- 
tions from the 1990 apparition, 108 observations from the 
1995 apparition, and 40 observations from the 2001 appari- 
tion. The composite lightcurve is shown in Fig. 11. This 
comet shows a post-perihelion excess but not so strong as 
C-G, d'Arrest, Tempel 2 or Borrelly. We can see an increas- 
ing scatter of the observations as the comet moves away from 
perihelion (the post-perihelion observations after about 60 
days were not considered for this reason). 

5.9 46P/Wirtanen 

This comet was the former target of the Rosetta mission, 
which will finally be encountering comet C-G in 2014. The 
lightcurve dataset for this comet contains 105 observations 
from the 1991 apparition. The lightcurve is shown in Fig. 12. 
We also studied the lightcurve corresponding to the 1997 ap- 
parition (with a dataset of 104 observations). We note that 
the 1997 observations seems to be somewhat brighter (about 
several tenths of magnitude) than those from 1991, so we de- 
cided not to combine both apparitions into one single com- 
posite lightcurve, since it is not clear that the lightcurve has 
not changed from one to the other. Since the 1997 observa- 
tions are much more scattered than those from 1991 as we 
go far from the perihelion, which results in an more poorly 
defined lightcurve, we decided to take the 1991 lightcurve 
as the nominal one. As inferred from this lightcurve, besides 
a post-perihelion brightness excess, this comet also shows a 
rapid decrease of the gaseous activity for increasing heliocen- 
tric distances. This is in agreement with Jorda and Rickman 
(1995), who also studied the 1991 lightcurve of this comet. 

5.10 67P/Churyumov-Gerasimenko 

This comet will be encountered by the Rosetta spacecraft 
in 2014, which will orbit it during its journey through the 
inner solar system, and land a probe on its surface. Hence 
a great wealth of physical data for this particular comet is 
expected in the next years. 

The lightcurve dataset for this comet contains 260 ob- 
servations from the 1982 apparition, and 127 observations 
from the 1996 apparition. The composite lightcurve is shown 
in Fig. 13. This comet shows a strong post-perihelion asym- 
metry. We also note that this comet is the faintest one of the 
sample (followed closely by Tempel 1), since its maximum 
brightness is somewhat below 10 mag, while most JF comets 
show maxima at magnitudes between 7 - 9, as we pointed 
out before. The post-perihelion excess and the maximum 
brightness obtained is in agreement with other works (e.g 
Ferrin 2005). 

5.11 81P/Wild 2 

This comet was visited by the Stardust spacecraft in January 
2, 2004, which obtained high resolution images of the nu- 
cleus, showing remarkable differences with previous images 
of Halley and Borrelly, like the evidence of impact craters 
(Brownlee et al. 2004). 

The lightcurve dataset for this comet contains 78 ob- 
servations from the 1990 apparition, 626 observations from 
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the 1997 apparition, and 72 observations from the 2003 ap- 
parition. The composite lightcurve is shown in Fig. 14. This 
comet shows the most strong pre-perihelion asymmetry of 
the sample. The bulk of the observations which defines the 
lightcurve around the maximum brightness corresponds to 
the 1997 apparition, while the 2003 observations contributes 
to define the pre-perihelion branch, and the 1990 observa- 
tions help to define the outermost part of the post-perihelion 
branch. This is the only comet of the sample which has 
a visual photometric coverage up to heliocentric distances 
of about 3 AU. Our lightcurve is in good agreement with 
that studied by Sekanina (2003). Ferrfn (2005) also derived 
a similar lightcurve for this comet, except that ours shows 
a steeper post-perihelion slope, which implies a more rapid 
fading of the activity. 



6 THE RESULTS 

For the computation of the non-gravitational change AP 
and the comet mass (from eqs. ([5]) and (0, respectively) 
we use some routines from Press et al. (1992) for numerical 
integration of functions and for root finding, and write our 
own routine to solve the Kepler equation. 

We first introduce the main sources of uncertainty in the 
computation of the cometary masses, and then the results 
obtained for the comets of the sample. Our intention was 
not to proceed to a rigorous computation of errors, since 
that would be meaningless under the general assumptions 
we made and the statistical scope of this work, but instead 
to evaluate in an empirical way the effect of the main sources 
of uncertainty in the computed masses. 

6.1 Sources of uncertainty 

The most important sources of uncertainty are the comet 
lightcurve, the model parameters (particularly the effec- 
tive outgassing velocity, which maybe is the most uncer- 
tain of the parameters assumed), and the computed non- 
gravitational change AP. 

6.1.1 The lightcurve 

In order to estimate the uncertainty due to the lightcurve, 
we proceeded as follows: as we explained in Section 4.1, the 
lightcurve is defined by a least-squared fitted polynomial to a 
set of selected observations. This set of observations was ob- 
tained by dividing the observational time interval in bins of 
fixed size, and selecting, in each bin, the three brightest ob- 
servations with the condition that they do not significantly 
depart from the bulk of observations. By varying the bin size 
(from 2 to 6 days, with one day step) , we picked different sets 
of brightest magnitudes leading to different lightcurves (one 
per bin size), and hence different cut-off magnitudes (rrihc)- 
For each one of these lightcurves (with its corresponding cut- 
off magnitude) we computed the mass (leaving the rest of 
the parameters involved fixed with their nominal values). In 
this way we obtained a set of different mass estimates for 
each comet of the sample. Hence, for a given comet, we take 
the mean value of the relative differences in the computed 
masses, with respect to the nominal value, as an empirical 



estimate of the uncertainty in the comet mass due to the 
lightcurve uncertainty. 

Comet d'Arrest was found to be the comet with the 
smaller uncertainty due to its lightcurve ((AM /M)n g f l tcurve 
= 3 %), followed by, in increasing order of uncertainty: Hal- 
ley and C-G (6 %), Tempel 1 and Tempel 2 (8 %), Kopff 
(9%), Encke (11 %), Borrelly (13 %), Wild 2 (28 %), H-M-P 
(43 %) and Wirtanen (44 %). 

We note that these values do not necessarily reflect the 
quality of the lightcurve in all cases, since, for instance, Tem- 
pel 1 has a relatively poorly defined post-perihelion branch 
of the lightcurve, though we got a small uncertainty in the 
computed mass. Nevertheless, in most cases, we find a good 
correlation between the lightcurve quality and the estimated 
uncertainty in the mass, as in the case of Wirtanen. This 
comet shows the most poorly defined lightcurve of the sam- 
ple, due to a relatively low number of observations, which 
reflects in the largest estimated uncertainty in the mass. 

6.1.2 The model parameters 

From the discussion of Section 3, we assumed an error 
of ± 0.1 km s _1 in < u >. We also assumed 0.05 ^ 
< | sin(r;) cos(</>)| > ^ 0.15, and 0.9 < < cos(?7) > < 1. 

Although a thorough study of the correlation between 
water production rates and visual magnitudes is beyond 
the scope of this work, we perform individual calibrations 
for some comets of our sample (as explained in Section 
4.2), just to evaluate, in an empirical way, how the uncer- 
tainty related to the constants {oi, 02} propagates to the 
computed masses. We found that the relative difference 
in the mass computed using the individual calibration, 
with respect to the nominal mass (i.e. the mass computed 
using the Jorda et al. (2008) calibration) varied between 
13 % (for Wild 2) and 46 % (for Tempel 1). The average 
of these relative differences was about 30 %. We also found 
that the relative difference in the mass computed using the 
Jorda et al. (1992) calibration, with respect to the nominal 
mass, was about 30 % for the comets of the sample, except 
for comet Halley for which the difference was about 20 
%. Hence, we estimate the uncertainty in the computed 
masses due to the calibration constants {01, 02} to be about 
(AM /M) calibration — 30 %, for all the comets of the sample. 



6.1.3 The non-gravitational effect 

For some comets we found in the literature an estimate of 
the absolute error A(AP) (e.g. Borrelly by Davidsson and 
Gutierrez 2004, Tempel 1 by Davidsson et al. 2007, C-G by 
Davidsson and Gutierrez 2005, and Wild 2 by Davidsson 
and Gutierrez 2006). For the rest of the comets we have to 
adopt some criteria in order to estimate the uncertainty of 
this parameter. For those comets which have showed some 
variation in the non-gravitational effect AP during the ap- 
paritions considered in this work (e.g. Encke, Kopff, H-M-P, 
Wirtanen), we estimate the error as the standard deviation 
with respect to the mean value (although for these comets 
intrinsic variations of the change AP cannot be ruled out, 
we neglect them by considering a composite lightcurve) . For 
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Figure 4. Heliocentric total visual magnitudes as a function of time (days from the perihelion passage), for 
comet lP/Halley. The polynomial fit m h (t) to the upper envelope of the broad distribution of photometric 
measurements is shown. The horizontal line indicate the cut- off magnitude m^c, as defined in Section 5. The 
value of rrihc is a l so shown. The magnitude data for this comet and the rest of the sample were extracted 
from Kamel 1991a (for observations before 1990) and from the ICQ archive (for observations from 1990 and 
forward) . 
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Figure 5. Idem as Fig. 4 for comet 2P/Encke. 
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Figure 6. Idem as Fig. 4 for comet 6P/d'Arrest. 
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Figure 7. Idem as Fig. 4 for comet 9P/Tempel 1. 



comets Halley and d" Arrest, in the absence of a better cri- 
terion, we assume an error given by the average of the nor- 
malized (i.e. relative) errors estimated for the former comets. 
For comet Tempel 2, which shows a very small AP like Tem- 
pel 1, we assumed the same error than for this latter comet. 



The estimated uncertainties range from ~ 4 % for comets 
Kopff and H-M-P, up to ~ 58 % for comet Wild 2. The 
intermediate values of the estimated uncertainty are: ~ 6 
% (comets Halley, d' Arrest and Borrelly), ~ 13 % (comet 
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Figure 8. Idem as Fig. 4 for comet lOP/Tempel 2. 
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Figure 9. Idem as Fig. 4 for comet 19P/Borrclly. 



Wirtanen), ~ 20 % (comet C-G), ~ 27 % (comets Tempel 
1 and Tempel 2), and ~ 50 % (comet Encke). 



6.1.4 Estimate of the overall error in the mass and 
density 

Although for some parameters involved in the computation 
of the mass (like the non-gravitational effect AP, and the 
effective outflow velocity < u >) their effect in the mass can 
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Figure 11. Idem as Fig. 4 for comet 45P/H-M-P. 



be easily evaluated separately, since the mass is proportional 
to them, we found more convenient to evaluate the overall 
uncertainty in the computed mass in an empirical way. In 
this regard we considered, for each comet of the sample, 
iV points of a parametric space defined by the parameters 



< u >, < cos(ry) >, < j sin(r?) cos(^>)| >, AP, and F, where 
F is a factor that accounts for the uncertainty due to the 
lightcurve and the calibration constants {01, 02}, defined as 
F = 1 + (AM/M), 

ightcurve 

+ (AM/M) calibration (a value of 
F = 1 implies a negligible uncertainty due to the lightcurve 
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and the calibration constants). We chose N = 1000. We 
denote by Fl the computed values of F for the comets of 
the sample. We obtained values ranging between Fl = 1.3 
for comet d'Arrest and Fl = 1.7 for comets H-M-P and 
Wirtanen. Each point of the parametric space was obtained 



by generating a random value with an uniform distribution 
for each parameter, with the constrains: 0.17 ^ < u > ^ 
0.37, 0.9 ^ < cos(t7) > < 1, 0.05 sC < | sm{rf) cos(0)| > 
«S 0.15, AP roi „ AP sC AP max , and 1 F «S F L , where 
AP min = AP nom -A(AP), and AP max = AP nom +A(AP), 
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according to the estimated absolute error A( AP) (computed 
from the relative errors presented in Section 6.1.3), and to 
the nominal value AP nom for the non-gravitational effect 
given in Table 3. Then, each one of these N points represents 
a set of parameter values from which a mass is computed, 
by means of eq. (J5]l , where the factor F multiplies the right 
hand member of this equation. By computing the standard 
deviation of these N mass estimates with respect to the 
nominal value we evaluate the overall uncertainty or error 
in the mass (AM) of a given comet. 

Finally, by propagating the estimated errors in the mass 
(AM) and in the effective nucleus radius (AP) (given in 
Sections 6.2.1 - 6.2.11), we obtain the error in the mass 
density 



A P 



4tt 



AM MAP 
+ 



P3 ' P4 

where P represents the effective nuclear radiufl 



(11) 



6.2.1 lP/Halley 

Mass. We derive a mass of [3.2 ±1.2] x 10 14 kg for this 
comet. 

Size and density. This comet has a determination of its nu- 
cleus size measured on close-up images from fly-by missions, 
which is the best procedure to determine unambiguously the 
size, albedo, or shape of a comet nucleus. Keller et al. (1987) 
determined an effective radius of 5.2 ± 1.0 km. According 
to this value for the comet size we derive a bulk density of 
0.5 ± 0.3 g cm -3 for this comet. 

Other works. Rickman (1989) estimated M ~ [1.3 - 3.1] x 
10 14 kg. Later, Skorov and Rickman (1999) revised the Rick- 
man's (1989) estimate and obtained p = [0.5 - 1.2] g cm -3 , 
with a preference for the lower value. Sagdeev et al. (1987) 
derived M = 3 x 10 14 kg, and a bulk density of 0.6±°;g g 
cm~ 3 , also based on a non-gravitational model. Peale (1989) 
estimated p = [0.03 - 4.9] g cm -3 , with a preference for p = 
0.7 g cm -3 . These results are consistent with ours. 



6.2 Computed masses and densities 

We present the computed masses from the nominal values for 
the model parameters (summarized in Table 1), with a rough 
estimate of the mass uncertainty due to the main sources of 
uncertainty as considered before. We also show the derived 
bulk mass density with an error estimate as derived from the 
size and the mass uncertainties, respectively. The results are 
compared with those from other works, and summarized in 
Table 5. 

2 The effective nuclear radius defines the radius of a sphere whose 
volume is equal to that of the comet nucleus (usually modeled as 
a triaxial ellipsoid of semiaxes a, b, c) 



6.2.2 2P/Encke 

Mass. We derive a mass of [9.2 ± 5.8] x 10 13 kg for this 
comet. 

Size and density. Encke has the shortest orbital period 
among the known comets, ft is also a NEC (i.e. a Near Earth 
Comet). As a consequence of its short periodicity, proximity 
to Earth, and orbital stability, this comet has been exten- 
sively studied from ground-based observations, and hence 
has several size determinations, also with different methods. 
Luu and Jewitt (1990) estimate an effective nucleus radius 
Pjve// of 3.28 ± 0.06 km from near-aphelion and time se- 
ries CCD photometry. Lowry and Weissman (2007) deter- 
mine R-Neff = 3.95 ± 0.06 km by using the same technique. 
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From radar and IR study Nolan and Harmon (2005) esti- 
mate R Ne ff = 2.42 to.H km. Tancredi et al. (2006) derive 
RNeff = 1-95 lg 67 km from the estimated absolute nuclear 
magnitude. Since there are for this comet several different 
estimates of the nucleus size, for which some of them yield 
unphysical values of p when the corresponding errors arc 
considered, we preferred to assume a sort of averaged radius, 
with a certain error estimate. After analyzing the different 
estimates, we adopted RNeff ~ 3 ± 1 km. The correspond- 
ing bulk density would be then p = 0.8 ± 0.8 g cm~ 3 . 
Other works. Rickman et al. (1987) estimated M ~ [2.4 - 
3.2] x 10 13 kg, which is about half the value we obtained. 
We find that the cause of this discrepancy could be that 
Rickman et al. considered a non-gravitational AP that was 
a factor of 2 larger than ours. 

6.2.3 6P/d' Arrest 

Mass. We derive a mass of [2.8 ± 0.8] x 10 12 kg for this 
comet. 

Size and density. Tancredi et al. (2006) estimate RNeff = 
1.7 ± 0.1 km. According to this value for the comet size we 
derive a bulk density of 0.15 ± 0.05 g cm -3 for this comet. 
Other works. Szutowicz and Rickman (2006) estimate M — 
7.0 x 10 12 kg, from non-gravitational modeling. This value 
is larger than ours by a factor of 2.7. The lightcurve ob- 
tained and the AP considered by these authors are similar 
to ours, so we looked for differences in the values assumed for 
the model parameters. Szutowicz and Rickman determine an 
empirical value for the lag angle n by fitting an asymmetric 
model to astrometric observations performed between 1851 
and 2001. They obtained n ~ 14°, which agrees very well 
with our assumption. On the other hand, they assume a 
mean effective outflow velocity < u > = 0.6 km s _1 , which 
is larger than the value we assume by a factor of 2.2. This 
would imply a very collimated flux with a very high mo- 
mentum transfer efficiency. Since the mass is directly propor- 
tional to this parameter, this explains the difference between 
both mass estimates. Szutowicz and Rickman also used dif- 
ferent calibration coefficients: they assumed the value given 
by Jorda et al. (1992) for the slope coefficient (ai = -0.240), 
and determine ai — 30.5 by fitting the observations to mea- 
sured gas production rates. We conclude that the different 
values assumed for the mean effective outflow velocity, as 
well as for the calibration coefficients, account for the dis- 
crepancy between both results. In principle, Szutowicz and 
Rickman's adopted value of < u > seems to be extremely 
high. 

6.2.4 9P/Tempel 1 

Mass. We derive a mass of [2.3 ±1.6] x 10 13 kg for this 
comet. 

Size and density. A'Hcarn et al. (2005) estimate RNeff = 
3.0 ±0.1 km from close-up images taken by the Deep Impact 
fly- by spacecraft. According to this value for the comet size 
we derive a bulk density of 0.2 ± 0.1 g cm~ 3 for this comet. 
Other works. Richardson and Melosh (2006) obtain M — 5.0 
x 10 13 kg by studying the expansion of the base of the coni- 
cal ejecta plume created by the impact of a 370 kg projectile 
delivered by the Deep Impact spacecraft on the comet nu- 
cleus. This is the most direct measurement obtained to date 



(and independent of the non-gravitational effect, as we would 
like to emphasize). They obtain a bulk density p = 0.4 ± 0.3 
g cm~ 3 . Davidsson et al. (2007) obtain M = (5.8 ± 1.6) x 
10 13 kg, and p — 0.45 ± 0.25 g cm~ 3 , by utilizing a sophisti- 
cated thermophysical model of the non-gravitational forces 
(which in addition relies on the measured gas production 
rates instead of the lightcurve). Our results are consistent 
with both independent works, which gives support to the 
usefulness of the modeling method used here. 

6.2.5 lOP/Tempel 2 

Mass. We derive a mass of [3.5 ±1.5] x 10 14 kg for this 
comet. 

Size and density. Jewitt and Luu (1989) obtain RNeff = 5.3 
lg' 2 km from time series CCD photometry. Tancredi et al. 
(2006) estimate i?jv e / / = 4.02 1^96 km from the estimated 
absolute nuclear magnitude. From the different estimates we 
adopted an averaged value RNeff ~ 4.8 ± 0.7 km. According 
to this value for the comet size, we derive a bulk density of 
0.7 ± 0.4 g cm~ 3 . 

Other works. Rickman et al. (1991) obtain M = (1.6 ± 0.5) 
x 10 14 kg, from non-gravitational force modeling. This is 
~ 50 % of the value derived in this work. The main cause 
of this discrepancy could be that Rickman et al. computed 
the mass considering only the radial non-gravitational force 
component, and hence neglecting the contribution from the 
transverse component (according to our results, the trans- 
verse component contributes to ~ 50 % of the mass esti- 
mate) .From the rotational lightcurve Jewitt and Luu (1989) 
derived parameters a/6 and P ro t (the axial ratio and the ro- 
tational period, respectively), and obtained a minimum bulk 
density p m in by assuming a negligible tensile strength of the 
nucleus material. They calculated pmin = 0.3 g cm~ 3 , which 
is consistent with our result. 



6.2.6 19P/Borrelly 

Mass. We derive a mass of [2.7 ± 2.1] x 10 12 kg for this 
comet. 

Size and density. From HST observations Lamy et al. (1998) 
obtain semi-axes a — A A km and 6 = c = 1.8 km by approxi- 
mating the nucleus to a triaxial ellipsoid, which corresponds 
to RNef f = 2.4 ± 0.2 km. This was later confirmed by direct 
imaging of the nucleus by Deep Space 1 (Soderblom et al. 
2002). According to this value for the comet size we derive 
a bulk density of 0.05 ± 0.04 g cm -3 for this comet. 
Other works. Davidsson and Gutierrez (2004) constrain the 
mass within the range M = [8 - 24] x 10 12 kg, and the 
bulk density in the range p = [0.10 - 0.30] g cm~ 3 , based 
on a thermophysical modeling of the non-gravitational force. 
Rickman et al. (1987) derive an upper limit M < 8.4 x 10 12 
kg, which is consistent with our result. Davidsson (2001) 
obtains a minimum bulk density p m in = 0.04 g cm -3 by re- 
quiring that the nucleus withstands disruption due to cen- 
trifugal forces, assuming a strengthless nucleus material. By 
comparison with these results, we may conclude that we 
somewhat underestimated the mass and bulk density of Bor- 
relly. Nevertheless, all the results suggest that we are deal- 
ing with an extremely low density comet. In disagreement 
with these results, Farnham and Cochran (2002) obtain a 
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mass estimate that is substantially higher: they derive M 
= 3.3 x 10 13 kg and p = 0.49 g cm" 3 by following the 
method introduced by Rickman et al. (1987), but using a 
different vaporization model for the water production. They 
observed a strong primary (sunward) jet that was aligned 
with the nucleus spin axis (according to their derived pole 
position) . They also assumed an relatively small active area 
(4 % of the nucleus surface area), close enough to the pole 
that it would be illuminated almost continuously during 
the rotation of the nucleus. The combination of these two 
features (spin orientation and location and size of the pri- 
mary source of activity) will produce a strong collimation 
of the jet, making the non-gravitational force much more 
efficient (e.g. the non-gravitational force would produce the 
same non-gravitational effect AP acting on a larger nucleus 
mass, than a less efficient non-gravitational force acting on 
a smaller nucleus mass). They assumed < u > = 0.33 km 
s" 1 . 

6.2.7 22P/Kopff 

Mass. We derive a mass of [5.3 ± 2.2] x 10 12 kg for this 
comet. 

Size and density. Tancredi et al. (2006) estimate Rneff = 
1.8 ± 0.2 km from the estimated absolute nuclear magni- 
tude. According to this value for the comet size we derive a 
bulk density of 0.2 ± 0.1 g cm -3 for this comet. 
Other works. Rickman et al. (1987) derive an upper limit M 
< 2.9 x 10 13 kg, which is consistent with our result. 

6.2.8 45P/Honda-Mkros-Pajdusakova 

Mass. We derive a mass of 1.9 x 10 1 kg for this comet, with 
an upper limit of 5.4 x 10 11 kg. The lower limit turns out 
to be negative which is unphysical. Therefore the computed 
mass of this comet is very uncertain and should be taken 
with extreme caution. The main reason for this uncertainty 
is that the radial and transverse terms I r , It that intervene 
in the computation of the mass (cf. eq. [5} turn out to be 
of similar absolute values but different signs (see Table 5 
below), so small errors in and It can lead to either positive 
or negative values in (/,. + It). 

Size and density. Tancredi et al. (2006) estimate RNeff = 
0.33 ± 0.07 km from the estimated absolute nuclear magni- 
tude. This would be the smallest comet of the sample, and 
one of the smallest known. According to this value for the 
comet size we derive a bulk density of 1.2 g cm -3 for this 
comet. But we must emphasize that the size estimate is very 
uncertain, according to the referred authors. Because of the 
large uncertainty (not only in the nuclear size but also in the 
computed mass) we decided to exclude this comet from the 
general discussion between masses and densities presented 
in Section 7. 

Other works. Rickman et al. (1987) derive M = [1.5 - 5.4] 
x 10 11 kg, which is consistent with our result. 

6.2.9 46P/Wirtanen 

Mass. We derive a mass of [3.3 ± 2.3] x 10 11 kg for this 
comet. 

Size and density. Tancredi et al. (2006) estimate RNeff = 



0.58 ± 0.03 km from the estimated absolute nuclear magni- 
tude. According to this value for the comet size we derive a 
bulk density of 0.4 ± 0.3 g cm -3 for this comet. 
Other works. We have not found another mass estimate for 
this comet. 



6.2.10 67P/Churyumov-Gerasimenko 

Mass. We derive a mass of [1.5 ± 0.6] x 10 13 kg for this 
comet. 

Size and density. From HST observations Lamy et al. (2006) 
obtain RNeff = 1-98 ± 0.02 km. According to this value for 
the comet size we derive a bulk density of 0.5 ± 0.2 g cm -3 
for this comet. 

Other works. Davidsson and Gutierrez (2005) obtain M = 
[0.35 - 2.1] x 10 13 kg, and the bulk density in the range p = 
[0.10 - 0.60] g cm -3 , based on a sophisticated thermophys- 
ical modeling of the non-gravitational force. Our result is 
consistent with their work. 



6.2.11 8 IP/ Wild 2 

Mass. We derive a mass of 8.1 x 10 12 kg for this comet, with 
an upper limit of 2.5 x 10 13 kg. The lower limit turns out 
to be negative which is unphysical. Therefore the computed 
mass of this comet is very uncertain and should be taken 
with extreme caution. The main reason for this uncertainty 
is the same as for comet 45P/H-M-P. 

Size and density. Brownlee et al. (2004) and Duxbury et 
al. (2004) found that the comet can be approximated by a 
triaxial ellipsoid with semi-axes {a,b,c} = {2.75,2.00,1.65} 
± 0.05 km, from Stardust images. This corresponds to 
R N eff = 2.09 ±o.06 km - Howington-Kraus et al. (2005) 
performed a refined image analysis obtaining {a, b, c} = 
{2.607,2.002,1.350} ± 0.001 km, which corresponds to 
RNeff = 1.917 km. According to this value for the comet 
size we derive a bulk density of 0.3 g cm -3 with an upper 
limit of 0.8 g cm -3 . 

Other works. Davidsson and Gutierrez (2006) determined 
an upper limit M ^ 2.3 x 10 13 kg, and p ^ 0.60 g cm -3 , 
based on a sophisticated thermophysical modeling of the 
non-gravitational force. Our result is consistent with their 
work. 

In Table 5 we summarize our results, namely the com- 
ponents I t , the computed masses and bulk densities, the 
latter based on the radii found in the literature, also included 
in the table. 

According to the results of Table 5, the range of 
cometary masses of our sample cover several orders of mag- 
nitude, ranging from ~ 3-4xl0 14 kg (Halley and Tempel 2) 
down to ~ 2-3 x 10 11 kg (H-M-P, Wirtanen), while the densi- 
ties are almost all of the same order (tenths of g cm -3 ), and 
below ~ 0.8 g cm -3 with a mean value « 0.4 g cm -3 (we ex- 
cluded comet H-M-P due to its larger uncertainty in the size, 
as compared to the rest of the comets). This result is consis- 
tent with the cut-off value ~ 0.6 g cm -3 for cometary den- 
sities suggested by Lowry and Weissman (2003), Snodgrass 
et al. (2006). In the case of Borrelly we obtain en extremely 
low density of 0.05 g cm -3 (very close to the lower limit 
imposed by the centrifugal forces, according to Davidsson 
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Table 5. Masses and densities. 



Comet 


I r 


h 


M 




Rcf . □ 


p 




(molec. AU _1 ) 


(molec. AU" 1 ) 


(kg) 


(km) 




(g cm -3 ) 


lP/Halley 


+1.66X10 36 


+1.33X10 36 


3.2X10 14 


5.2 


(a) 


0.5 


2P/Encke 


-8.34X10 34 


-6.96X10 34 


9.2X10 13 


3 


(b) 


0.8 


6P/d' Arrest 


+2.56X10 34 


+1.09X10 34 


2.8X10 12 


1.7 


(c) 


0.15 


9P/Tempel 1 


-1.60X10 33 


+5.42 xlO 33 


2.3X10 13 


3.0 


(d) 


0.2 


lOP/Tempel 2 


+3.31 xlO 34 


+3.16X10 34 


3.5X10 14 


4.8 


(e) 


0.7 


19P/Borrelly 


+1.39X10 34 


-2.98X10 34 


2.7X10 12 


2.4 


(f) 


0.05 


22P/KopfT 


-1.55X10 34 


-3.24X10 34 


5.3X10 12 


1.8 


(c) 


0.2 


45P/H-M-P 


+2.45 xlO 34 


-3.26X10 34 


(1.9 xlO 11 ) 


0.33 


(c) 


(1.2) 


46P /Wirtanen 


+8.26X10 32 


-1.24X10 34 


3.3x10" 


0.58 


(c) 


0.4 


67P/C-G 


+1.46X10 34 


+6.40X10 33 


1.5xl0 13 


1.98 


(g) 


0.5 


81P/Wild 2 


-1.67X10 34 


+2.28X10 34 


(8.1xl0 12 ) 


1.917 


(t) 


(0.3) 



a (a) Keller et al. (1987), (b) See Section 6.2.2, (c) Tancredi et al. (2006), (d) A'Hearn et al. (2005), (e) See 
Section 6.2.5, (f) Lamy et al. (1998), (g) Lamy et al. (2006), (h) Howington-Krauss et al (2005). 



2001) though the large uncertainty of this estimate admits 
values more in line with the rest of the comets. Thermophys- 
ical modeling provides a higher value for Borrelly's density, 
though still pointing to a very low-density comet. 



6.3 The correlation between the non-gravitational 
change AP and the perihelion asymmetries of 
the lightcurves 

We can quantitatively evaluate the perihelion asymmetry 
of the gas production curves by introducing the following 
definition 

A = Q- 1 f Q(t)sm[f(t)]dt (12) 
Jo 

where t, f and P are defined as before, and Q m is a normal- 
ization factor arbitrarily chosen to be the maximum water 
production rate of the comet. According to this definition, A 
is measured in days. In Fig. 15 we plot the perihelion asym- 
metry against the time of the maximum observed brightness 
as derived from the lightcurves (we remind the reader that 
the time is referred to the time of perihelion passage) . As it 
should be expected, we can observe in general a linear re- 
lationship between both parameters, which also show very 
similar values. 

We find that for those comets which show a strong pos- 
itive perihelion asymmetry (i.e. a notorious post-perihelion 
excess of gaseous activity), the radial I r contribution domi- 
nates, with the exception of Borrelly. These are the cases of 
comets Halley, d'Arrest, Tempel 2, and C-G (actually, in the 
case of Tempel 2 the radial dominance is very weak). Wirta- 
nen shows the most symmetrical lightcurve of the sample, 
with a weak positive perihelion asymmetry. In this case the 
transverse term I t turns out to be the dominant one, as ex- 
pected. Comets Kopff and Wild 2 show a strong negative 
perihelion asymmetry (i.e. a notorious pre-perihelion excess 
of activity), and a dominant transverse contribution. Even 
though the strong asymmetric lightcurves cause large radial 
components in the latter two comets, their transverse contri- 
butions are actually large too (a factor ~ 1.4-2 greater than 
the radial ones) . Tempel 1 also shows a negative asymmetry 



though a very weak one, so the transverse term dominates. 
Encke shows an intermediate negative perihelion asymme- 
try, but a slight dominant radial contribution instead. 

Festou et al. (1990) and Rickman et al. (1991b) stud- 
ied the correlation between the perihelion asymmetries of 
the gas production curves of periodic comets and the non- 
gravitational perturbations of their orbital periods. In this 
regard they defined the standard non- gravitational change 
AP' as the non-gravitational effect in the orbital period that 
the comet would experience if its semimajor axis was 3.5 AU, 
which is given by 

AP' = AP(^)- (13) 

where a is the semimajor axis. AP' then represents the 
non-gravitational effect AP corrected for different semi- 
major axes. Their result was that, in general, the perihe- 
lion asymmetries give the dominant contribution to the non- 
gravitational effect of the orbital period. Actually they found 
a linear correlation between the perihelion asymmetry and 
AP'. This is expected if the I r contribution dominates, but 
not otherwise. 

In Fig. 15 we plot the perihelion asymmetry (as defined 
in eq. (|12[) , which is not exactly the same definition given by 
Festou et al. 1990, but it should be equivalent) as a function 
of AP' for the comets of the sample. According to Festou et 
al. (1990) and Rickman et al. (1991b), it should be expected 
an increase of the perihelion asymmetry with the increasing 
non-gravitational AP' (in absolute value). As it can be seen 
in the plot (bottom panel), we did not find such a linear 
correlation. As shown there, for some comets the transverse 
component It becomes dominant, which weakens or com- 
pletely blurs any possible correlation between AP' and the 
radial component I r . 



7 DENSITIES VS. MASSES IN THE SOLAR 
SYSTEM 

We analyze the results obtained for our comet sample plac- 
ing them into a broader context of different populations of 
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Figure 15. Top: Perihelion asymmetry A as function of the time of the brightness maximum T for the 
comets of the sample. A lineal fit to the computed values is also shown. Bottom: Perihelion asymmetry A 
as a function of the standard non-gravitational effect AP' for the comets of the sample. Those comets with 
a dominant contribution from the radial force component to the non-gravitational effect are marked with 
open triangles, and those comets with a dominant transverse force component contribution are marked with 
filled triangles. 



bodies in the Solar System. In Fig. 16 we plot the bulk den- 
sity as a function of the mass for the different populations 
of minor bodies in the Solar System as indicated in the fig- 
ure caption. In the case of the comets we plot the computed 
masses and densities, indicating the estimated errors for the 
computed masses and bulk densities as derived in Section 
6.2. Fig. 16 shows that the computed comets form a more 
compact group than the other populations of minor bodies, 
and also present the lowest densities (p < 0.8 g cm -3 ). 
Our computed comets and the NEAs occupy the left of the 
diagram, corresponding to the smaller masses (except one 
NEA). The main distinction between both populations in 
this space of parameters is that the NEAs of the sample have 
larger bulk densities. Only two NEAs show p < 1 g cm -3 , 
but still above the upper limit of ~ 0.8 g cm -3 derived for 
the comets. Of course due to the smallness of the samples 
and the large uncertainties involved in the computation of 
masses and densities, both populations could actually some- 
what overlap. This would be consistent with the idea that 
some NEAs could be deactivated JF comets (either extinct 
or dormant). Nevertheless, the clear separation in the den- 
sity domain shown by comets and NEAs suggests that most 
NEAs are bona fide asteroids in agreement with Fernandez 
et al. (2002). 

We find other objects not cataloged as comets but also 
with p < 1 g cm" 3 ; these are: Patroclus, Amalthea, 1999 
TC36, Antiope and Emma (the last two are MB As). But 
in these cases the masses are at least about four orders of 
magnitude greater than the computed comet masses. 



We note in Fig. 16 that no bodies with masses M it 
10 20 kg have bulk densities i$ 1.5 g cm -3 . If we assume that 
comets and other ice-rich bodies of the outer solar system 
have similar compositions consisting of dust-ice mixtures, 
then the different densities would suggest that compaction 
of the material by self-gravity has taken place, so fluffy ob- 
jects should not be expected for masses bigger than about 
10 20 kg. We also have in the figure rocky bodies (NEAs 
and main-belt asteroids) whose differences in density may 
be due in this case to the different chemical composition, 
though different degrees of porosity might also be allowed. 
This is the case of most main-belt asteroids, TNOs, and, of 
course, the Moon and the largest Jovian satellites. There- 
fore, self-gravity would lead to noticeable variations in the 
primordial physical structure, ranging from loose aggregates 
to compacted and hardened material, and hence to a wide 
range of densities. In the case of the studied comets, due 
to the ranges of masses and sizes involved, we do no ex- 
pect them to increase their densities by self gravity. Hence 
comets would be fluffy, low-density objects which have pre- 
served their primordial fragile and porous structure. 



8 CONCLUDING REMARKS 

From non-gravitational force modeling we derive masses 
and densities for ten short-period comets of known 
sizes: lP/Halley, 2P/Encke, 6P/d'Arrest, 9P/Tempel 1, 
lOP/Tempel 2, 19P/Borrelly, 22P/Kopff, 46P/Wirtanen, 
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Figure 16. Bulk density as a function of the mass for the different populations of minor bodies of the 
Solar System. For clarity of the figure, we only give error bars for our sample comets. The represented 
objects are: the ten comets of our sample (comet 45P/H-M-P is omitted due to its large uncertainty, and for 
comet 81P/Wild 2 only the nominal values and upper limits are given), nineteen objects from the Main Belt 
Asteroids (MBAs) population, six objects from the Near-Earth Asteroids (NEAs) population (Eros, Itokawa, 
1999 KW 4 , 2000 DP107, 2000 UG11, 2002 CE26), one Trojan of Jupiter (Patroclus), one Centaur (Ceto- 
Phorcys), three trans-neptunian objects or TNOs (Pluto/Charon, 1999 TC36, 2003 ELgi), three terrestrial 
satellites (Moon, Phobos and Deimos), and ten Jovian satellites. The mass and density (or size) data for 
these objects comes from: Yeomans et al. (2000) (433 Eros), Shinsuke et al. (2006) (Itokawa), Merline et al. 

(2002) (1999 KW4, 2000 DP107, 2000 UGn), Shepard et al. (2006) (2002 CE 26 ), Kovacevic and Kuzmanoski 
(2007) (Ceres), Thomas et al. (2005) (size of Ceres and Pallas), Goffin (2001) (Pallas), Viateau and Rapaport 
(2001) (Vesta, Phartenope), Michalak (2001) (Hygiea, Eunomia), Lupishko (2006) (Psyche), Marchis et al. 

(2003) (Kalliope), Marchis et al. (2005b) (Eugenia, Camilla, Elektra, Emma, Huenna), Marchis et al. (2005a) 
(Sylvia), Marchis et al. (2005c) (Hermione), Descamps et al. (2005) (Antiope), Merline et al. (2002) (Ida, 
Pulcova), Marchis et al. (2006) (Patroclus), Grundy et al (2007) (Ceto-Phorcys), Tholen and Buie (1997) 
(Pluto/Charon), Stansberry et al. (2006) (1999 TC 36 ), Rabinowitz et al. (2006) (2003 EL 6 i), Anderson et 
al. (2005) (Amalthea), and references in De Pater and Lissauer (2001) (Moon, Phobos, Deimos, Amalthea, 
Io, Europe, Ganymedes, Calysto, Titan, Titania, Oberon, Triton, Nereida). 



67P/Churyumov-Gerasimenko and 81P/Wild 2. For another 
comet, 45P/H-M-P, it was not possible to compute a reliable 
mass and density. Our procedure follows the pioneer work of 
Rickman and colleagues (e.g. Rickman 1986, 1989; Rickman 
et al. 1987), and confirms the importance of the knowledge 
of non-gravitational effects on the comet's motion, as a fun- 
damental tool to derive cometary masses. The present work 
is based on lightcurve data, the non-gravitational term A2, 
and different assumptions for some physical parameters. Our 
main conclusions can be outlined as follows: 

(i) The derived cometary masses cover more than three 
orders of magnitude, within the range ~ [0.3 - 400] x 10 12 
kg. Our best estimate for the bulk densities are below 0.8 
g cm -3 , which is in agreement with models of the comet 
nucleus depicting it as a very low density object. The mean 
bulk density of our sample is 0.4 g cm -3 . 

(ii) The results obtained are in general consistent with 
other works. Particularly, we remark the consistency with 



results from sophisticated thermophysical non-gravitational 
force modeling, and from the Deep Impact team using a 
different method for Tempel 1. In particular, this agree- 
ment gives us some reassurance of the usefulness of the 
non-gravitational method used here for deriving cometary 
masses, despite all the uncertainties inherent to it. 

(iii) The method also relies on the empirical calibra- 
tion between visual magnitudes and water production rates, 
which, in spite of its poor physical understanding, it has 
proved to be a very useful tool. Since there is an ample and 
comprehensive database of total magnitudes, the method 
could be extended to many more periodic comets, for which 
the non-gravitational change AP has been computed. If 
we also know comet sizes (whose knowledge is also rapidly 
improving) we can derive the comet's bulk density. Even 
though the computed bulk density for a given comet has a 
large uncertainty, by computing them for a large sample of 
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comets we could narrow down the uncertainty of the whole 
sample to meaningful values. 

(iv) In connection with the previous points, we want 
to stress that our results are meaningful from a statistical 
sense, by considering an average value of a large sample of 
comets. Of course, this later statement is true if there are 
not systematic errors in the method by overestimating or 
underestimating any of the parameters used in the compu- 
tation of the masses. Even tough this possibility cannot be 
ruled out, we have been very careful in choosing the values 
of the input parameters within physically realistic ranges. 
We should warn, however, that, given the large errors bars, 
we cannot tell at this point whether there is a real dispersion 
in the bulk densities of comets reflecting different degrees of 
compaction of the material and/or differences in the chemi- 
cal composition. 

(v) We conclude that a strong perihelion asymmetry of 
the lightcurve should not be necessarily related to a negli- 
gible transverse force component contribution to the non- 
gravitational change in the orbital period. Actually, in some 
cases (like comets Kopff and Wild 2), the transverse com- 
ponent is found to be more important than the radial con- 
tribution. 

(vi) Limitations of the method. The mass solutions from 
eq. ([SJl are quite robust to small changes in the input pa- 
rameters except when the terms I r and It turn out to be of 
similar absolute value but different signs (as in the cases of 
45P/H-M-P and 81P/Wild 2). 

(vii) Sources of uncertainty. The lightcurve uncertainty 
and the values assumed for the model parameters, particu- 
larly for the effective outflow velocity < u >, would be in 
general the most important contribution to the uncertainty 
in the computed masses, although in some cases the non- 
gravitational change AP could become a major source of 
uncertainty too. An underestimate or overestimate of, e.g. 
< it >, would lead to a systematic underestimate or overes- 
timate of the computed mass. A more detailed study of this 
parameter, like its dependance on the heliocentric distance, 
has to wait for more elaborated results from thermophysical 
models as well as more observations, particularly in situ. 
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